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We study the strong coupling limit of the 2-flavor massless Schwinger model on a lattice using 
staggered fermions and the Hamiltonian approach to lattice gauge theories. Using the correspon- 
dence between the low-lying states of the 2-flavor strongly coupled lattice Schwinger model and the 
antiferromagnetic Heisenberg chain established in a previous paper, we explicitly compute the mass 
spectrum of this lattice gauge model: we identify the low-lying excitations of the Schwinger model 
with those of the Heisenberg model and compute the mass gaps of other excitations in terms of vac- 
uum expectation values (v.e.v.'s) of powers of the Heisenberg Hamiltonian and spin-spin correlation 
functions. We find a satisfactory agreement with the results of the continuum theory already at the 
second order in the strong coupling expansion. We show that the pattern of symmetry breaking 
of the continuum theory is well reproduced by the lattice theory; we sec indeed that in the lattice 
theory the isoscalar (Mi'ip) an( i isovector ($a a il)) chiral condensates are zero to every order in the 

strong coupling expansion. In addition, we find that the chiral condensate < ip^ ^r^r > 
is non zero also on the lattice; this is the relic in this lattice model of the axial anomaly in the 
continuum theory. 

We compute the v.e.v.'s of the spin-spin correlators of the Heisenberg model which are pertinent 
to the calculation of the mass spectrum and obtain an explicit construction of the lowest lying states 
for finite size Heisenberg antiferromagnetic chains. 



1 Introduction 



One of the analytical approaches to the study of gauge theories with confining spectra is the 
strong coupling expansion. In the strong coupling limit, confinement is explicit, the confining 
string is a stable object H and some other qualitative features of the spectrum are easily obtained. 
Formulation of the strong coupling expansion requires a gauge invariant ultraviolet cutoff, which is 
most conveniently implemented using a lattice regularization. 

One of the most difficult problems of the strong coupling approach to lattice gauge theory is its 
extrapolation to the continuum limit, which usually occurs at weak coupling. One symptom of this 
problem is that many choices of strong coupling theory produce identical continuum physics. In 
spite of this difficulty, there are strong coupling computations which claim some degree of success g, 
||, ||, ||, H . It is well known that between strongly coupled lattice gauge theories and quantum spin 
systems there is an intimate relationship |?], |^, ^) which is most useful to analyze chiral symmetry 
breaking in a variety of lattice models. 

A useful test of strong coupling expansions in lattice gauge theories is the study of models whose 
solution in the continuum is known even in the strong coupling regime. For example, the 1-flavor 
Schwinger model has been studied using the Hamiltonian lattice field theory within the strong 
coupling expansion ||, ^j; several relevant parameters, analytically computed on the lattice, are 
shown to be in good agreement with those of the exact continuum solution Q]. Furthermore, issues 
involving the realization of chiral symmetry on the lattice can be studied systematically . The 

1- flavor Schwinger model has also been used as an example of the fact that quantum link models 
may reproduce the physics of conventional Hamiltonian lattice gauge theories ]lf| . The physically 
more interesting case of the 2-flavor Schwinger model has been studied in [pTpl; in addition to 
the issues of spectrum and chiral symmetry breaking, one sees that the strong coupling limit of 
this model is mapped |H], | on a relevant quantum spin model — the one-dimensional spin-1/2 
quantum Heisenberg antiferromagnet. This correspondence is very useful since the ground state 
of the antiferromagnetic chain is known |l2] ] and its energy has been computed in M]; moreover, 
the complete spectrum has been determined by Faddeev and Takhtadzhyan using the algebraic 
Bethe ansatz. 

There are by now many hints at a correspondence between quantized gauge theories and quantum 
spin models, aimed at analyzing new phases relevant for condensed matter systems fl5| . Recently, 
Laughlin has argued that there is an analogy between the spectral data of gauge theories and 
strongly correlated electron systems Moreover, certain spin ladders have been related to the 

2- flavor Schwinger model in Ref. jll] where a relation between the physical parameters of the spin 
and the gauge systems was also found. 

In this context, the mapping between the strongly coupled 2-flavor Schwinger model and the 
quantum Heisenberg antiferromagnetic chain provides a concrete computational scheme in which 
the issue of the correspondence between quantized gauge theories and quantum spin models may 
be investigated. Because of the dimensionality of the coupling constant in (l+l)-dimensions, the 
infrared behavior is governed by the strong coupling limit, and it is tempting to conjecture the 
existence of an exact correspondence between the infrared limits of the Heisenberg and 2-flavor 
Schwinger models. 

In this paper we revisit the strong coupling limit of the 2-flavor lattice Schwinger model in the 
Hamiltonian formalism with staggered fermions. Using the results of Ref. Q we analyze in detail the 
spectrum of the model: the gauge theory vacuum is the ground state of the Heisenberg antiferromag- 
netic chain and we find that the low-lying excitations of the gauge model have the same quantum 
numbers of the states of the spin chain. In addition we compute the masses of the excitations to 
the second order in the strong coupling expansion; as pointed out in the computation needs the 
knowledge of some spin-spin correlators of the quantum Heisenberg antiferromagnetic chain, which 
we explicitly compute (l7], [l8) . Our analysis hints to the existence of two other massive isotriplet 
states in addition to the expected pseudoscalar massive isosinglet. 

We analyze then the pattern of symmetry breaking in the lattice theory. Even though the 
continuum axial symmetry is broken explicitly by the staggered fermions, a discrete axial symmetry 
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remains, it corresponds to a chiral rotation of 7r/2 and appears in the lattice theory as a translation 
by one lattice site. Since the ground state of the Heisenberg chain is unique and translationally 
invariant - at variance with the one flavor Schwinger model - in the 2-flavor model the discrete 
axial symmetry is unbroken. The chiral anomaly in this model is realized on the lattice via the 
explicit breaking of the J7 J 4(l)-symmetry induced by the staggered fermions. Therefore, the pattern 
of symmetry breaking on the lattice reproduces faithfully the one of the continuum theory. 

We also compute the chiral condensates in the strongly coupled 2-flavor lattice Schwinger model 
showing that the results of the lattice are in agreement with those of the continuum. In the 
continuum there is neither an isoscalar (rjjtp) nor an isovector {^)a a ^p') chiral condensate, since this 
is forbidden by the Coleman theorem |nj . We show that on the lattice these fermion condensates are 
also zero to all the orders in the strong coupling expansion. We find that — just as in the continuum 

theory — the non- vanishing chiral condensate is given by the v.e.v. < ip L ^tp L ^ ^r^r >j which we 
compute up to the second order in the strong coupling expansion. 

We finally use both the mass spectrum and the non-vanishing chiral condensate to compute, 
by means of suitable "Pade approximants" , the physical parameters of this lattice model and then 
compare their numerical values with the exact results of the continuum theory. 

In section 2 we review known results |2l| of the continuum 2-flavor Schwinger model and 
define its Hamiltonian lattice version. We also identify the lattice counterparts of the relevant 
symmetries used in [^0| for the description of the spectrum of the continuum theory. 

In section 3 we expand the results of Ref. || and set up the formalism needed for the strong 
coupling analysis of the lattice 2-flavor Schwinger model. We show that the ground and excited 
states of the antiferromagnetic Heisenberg model have the spectrum and quantum numbers which 
are expected for the ground state and massless excitations of the 2-flavor Schwinger model. 

In section 4 we construct the operators that create the massive excitations when acting on the 
ground state. We set up the strong coupling expansion and compute the corrections to the energies 
of the ground state and of the low- lying massive bosonic excitations. Subtracting the energy of the 
ground state from those of the excitations defines the lattice meson masses. 

In section 5 we show that both the isoscalar and the isovector chiral condensates are zero to every 
perturbative order in the strong coupling expansion due to the translational invariance of the strong 
coupling ground state. We also compute, up to the second order in the strong coupling expansion, 

the chiral condensate < ipj^ipL ^r^r >) which is the order parameter for the breaking of the 
Ua{^) symmetry. 

Section 6 is devoted to the comparison of the lattice results with the continuum theory; there 
we shall see that the lattice theory, properly extrapolated to the continuum via the use of Pade 
approximants, well reproduces the parameters of the continuum theory already at the second order 
of the strong coupling expansion. 

Section 7 is devoted to some concluding remarks. 

In the appendix A, after reviewing known results about the Bethe Ansatz solution |Q of the 
antiferromagnetic Heisenberg chain, we study explicitly the complete spectrum of finite size Heisen- 
berg antiferromagnetic chains of 4 and 6 sites in order to compare it to the Bethe Ansatz solution 
obtained in the thermodynamic limit. We write down explicitly the ground states of 4,6 and 8 site 
chains and we find the that, even for these very small systems, the results exhibited in jlj] are very 
well reproduced. This provides also an useful intuitive picture of the strong coupling ground state 
of the 2-flavor lattice Schwinger model. 

In appendix B we comment on the computation of the spin-spin correlator of the Heisenberg 
chain and provide a link between the results given in w\ and the method used in |lq| . 
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2 Two flavor Schwinger model in the continuum and on the 
lattice 



The action of the 1 + 1-dimensional electrodynamics with two charged Dirac spinor fields is 

2 



d 2 x 



,a=l 



— F F^ v 

4e? ^ 



(2.1) 



The theory has an internal SUl{2) (g) SUr (2)-flavor isospin symmetry; the Dirac fields are an 
isodoublet whereas the electromagnetic field is an isosinglet. It is well known that in 1 + 1 dimensions 
there is no spontaneous breakdown of continuous internal symmetries, unless there are anomalies 
or the Higgs phenomenon occurs. Neither mechanism is possible in the 2-flavor Schwinger model 
for the SUl(2) S l J7fl(2)-symmetry: isovector currents do not develop anomalies and there are no 
gauge fields coupled to the isospin currents. The particles belong then to isospin multiplets. For 
what concerns the U(l) gauge symmetry there is an Higgs phenomenon [22]] . 

The action is invariant under the symmetry 



SU L {2) ® SU R (2) <g> U v (l) <8 U A (l) 
The group generators act on the fermion isodoublet to give 



SU L {2) 


lpa(x) - 


-> (e l9a ^ PL ) ab i>b{x) , Mx) 


^Tb (x)(e- i6 «^ p *) ba 


SU R (2) 


1pa(x) - 


-> {e w ^ PR ) ab M%) , Mx) 


— (e~ w ^ p ^) ba 


Uy(l) 


4>a(x) " 


- (e^ x)1 ) ab Mx) , 4>l{x) - 


-^(z) (e- iS ^) ha 


Mi) 


lpa(x) - 


- (e^U M*) > ^) ~ 





where a a are the Pauli matrices, 6 a , 9{x) and a are real coefficients and 



(2.2) 

(2.3) 
(2.4) 
(2.5) 
(2.6) 
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1(1 + 76) 



(2.7) 



At the classical level the symmetries (2^5— 2J3) lead to conservation laws for the isovector, vector 
and axial currents 



f a {x)R 
jg(x) 



^aix^PL^abMx) 

^a( x h >11 ab^ b (x) 
i>a( X h >J 'l 5l a.b1pb(x) 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



It is well known that at the quantum level the vector and axial currents cannot be simultaneously 
conserved, due to the anomaly phenomenon [^3|. If the regularization is gauge invariant, so that 
the vector current is conserved, then the axial current acquires the anomaly which breaks the 
V A ( 1 )-symmetry 



0/V> - 



9^35 (x) 



,F» v (x) 



(2.12) 



J 2^ 

The isoscalar and isovector chiral condensates are zero due to the Coleman theorem fTs} ]; in fact, 
they would break not only the Ua{\) symmetry of the action, but also the continuum internal 
symmetry SU L (2) ® SU R (2) down to SU V (2). There is, however, a SU L {2) ® SU R (2) invariant 
operator, which is non-invariant under the LOi(l)-symmetry; it can acquire a non-vanishing v.e.v 
without violating Coleman's theorem and consequently may be regarded as a good order parameter 



for the J7 J 4(l)-breaking. Its expectation value is given by |24, 25 



<f>^MVM ) >= 



'47T 7T C 



(2.13) 
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It describes a process in which two right movers are anihilated and two left movers are created. 
Note that F, being quadrilincar in the fields, is actually invariant under chiral rotations of 7r/2, 
namely under the discrete axial symmetry 



(2.14) 



As a consequence, this part of the chiral symmetry group is not broken by the non- vanishing v.e.v. 
of F (HI). 

As we shall see in section 6, the lattice theory faithfully reproduces the pattern of symmetry 
breaking of the continuum theory; this happens even if on the lattice the SU (2)-flavor symmetry is 
not protected by the Coleman theorem. The isoscalar and isovector chiral condensates are zero also 
on the lattice, whereas the operator F acquires a non-vanishing v.e.v. due to the coupling of left 
and right movers induced by the gauge field. The continuous axial symmetry is broken explicitly 
by the staggered fermion, but the discrete axial symmetry (2.14) remains. 



The action (2.1) may be presented in usual abelian bosonized form |20|. Setting 



1 



V" 

the electric charge density and the action read 

jo 



—e^d^ a , a = 1,2 



-t=0x(*i 

\/7T 



S = J d 2 x 
By changing the variables to 



2lT 



<J> 2 



(*l + * 2 ) s 



•*a) 



one has 



S = 



d 2 x 



S. 

7T 



(2.15) 

(2.16) 
(2.17) 

(2.18) 
(2.19) 

(2.20) 



The theory describes two scalar fields, one massive and one massless. $+ is an isosinglet as 
evidenced from Eq.( 2.16| ); its mass ms — \J^ e c comes from the anomaly Eq.(2.12) |22 . Local 
electric charge conservation is spontaneously broken, but no Goldstone boson appears because the 
Goldstone mode may be gauged away. represents an isotriplet; it has rather involved nonlinear 
transformation properties under a general isospin transformation. All three isospin currents can be 
written in terms of but only the third component has a simple representation in terms of $_; 
namely 



^i(i)7^(i)-^(x)7^(i) := (2tt) 



(2.21) 



The other two isospin currents j^(x) and j 2 (x) are complicated nonlinear and nonlocal functions of 
$_ p0| ]; a more symmetrical treatment of the bosonized form of the isotriplet currents is available 
within the framework of nobn abelian bosonization [ p6| . For the multiflavor Schwinger model this 
approach has been carried out in p7| , providing results in agreement with [ p0| . 

The excitations are most conveniently classified in terms of the quantum numbers of P-parity 
and G-parity; G-parity is related to the charge conjugation G by 



G 



^G 



(2.22) 
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$_ is a G-even pseudoscalar, while <&+ is a G-odd pseudoscalar 

: 7 PG = 1-+ (2.23) 
$+ : 7 PG = 0"" . (2.24) 



The massive meson $ + is stable by G conservation since the action (2.20) is invariant under <I> + — ► 

In the massive SU(2) Schwinger model — when the mass of the fermion m is small compared 
to e 2 (strong coupling) — Coleman ||cj showed that — in addition to the triplet $_ (I PG = 1~ + ) 
— the low-energy spectrum exhibits a singlet I PG = ++ lying on top of the triplet <&_. In this 
limit the gauge theory is mapped to a sine-Gordon model and the low-lying excitations are soliton- 
antisoliton states. When m — » 0, these soliton-antisoliton states become massless [Slj; in this limit, 
the analysis of the many body wave functions, carried out in Ref. (23], hints to the existence of a 
whole class of massless states with positive G-parity; P-parity however cannot be determined with 
the procedure developed in These are not the only excitations of the model: way up in mass 
there is the isosinglet I PG = , ($+); already discussed in Ref. The model exhibits also 

triplets, whose mass — of order ms or greater — stays finite among the triplets there is a 
G-even state Q 

The Hamiltonian, gauge constraint and non-vanishing (anti-) commutators of the continuum 
2-flavor Schwinger model are 



H= dx 



~^E 2 (x) + ELi ( id * + eA W) M*)] ( 2 - 2 5) 

d x E(x) + ELi V'Kx)Va(x) ~ (2.26) 

[^(a;) J £;(tf)]=i*(x-i/) ,{i(i),4(!/)}=W(i-!/) ■ (2.27) 

A lattice Hamiltonian, constraint and (anti-) commutators reducing to ( [2.25[ 2~26| ,2.27) in the naive 
continuum limit are 

2 N N 2 

x — 1 sc=l a— 1 

E x -E x - X + ^L^+^L^-l^O , (2.28) 



The fermion fields are defined on the sites, x = 1, ...,N, the gauge and electric fields, A x and -E^, 
on the links [x; x + 1], A is an even integer and, when TV is finite it is convenient to impose periodic 
boundary conditions. When A is finite, the continuum limit is the 2-flavor Schwinger model on a 
circle |2^]. The coefficient t of the hopping term in (2.2S) plays the role of the lattice light speed. 
In the naive continuum limit, ej, = e c and t = 1. 

The Hamiltonian and gauge constraint exhibit the discrete symmetries 

• Parity P: 

A x — » -A_ x _i, £ a — > Va,* — (-l)^a,-., — » (-l)"V4,-x (2-29) 

• Discrete axial symmetry Y: 

A x ► A x+1 , E x > B^+x, Va,^ — > ip a ,x+i, i>l, x — ► ^i, x +i (2-30) 

• Charge conjugation C: 

A x — > -A x+ i, E x — > -E x+1 , if)a }X — > ipl x+1 , ipLx — > $a,x+i (2.31) 



X K. Harada private communication. 
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• G-parity: 



A x — ► — Ax+i, E x — > —E x+ \ 
i>i,x — > $2,x+l> V4,x — ► 4>2,x+l 



(2.32) 



The lattice 2-flavor Schwinger model is equivalent to a one dimensional quantum Coulomb gas 
on the lattice with two kinds of particles. To see this one can fix the gauge, A x = A (Coulomb 
gauge). Eliminating the non-constant electric field and using the gauge constraint, one obtains the 
effective Hamiltonian 



Hs — H u + H„ = 



■r,U 



+ 



it 
2a 



where the charge density is 
and the potential 



JV-1 



a i2-Kn(x—y) / N 



4 sin 2 ^ 



(2.33) 

(2.34) 
(2.35) 



is the Fourier transform of the inverse laplacian on the lattice for non zero momentum. The constant 
modes of the gauge field decouple in the thermodynamic limit N — ► oo. 



3 The strong coupling limit and the antiferromagnetic Heisen- 
berg Hamiltonian 

In a previous paper |(| we showed that the low-lying spectrum of the 2-flavor lattice Schwinger 
model in the strong coupling limit is equivalent to the spectrum of the antiferromagnetic Heisenberg 
model. There we also showed that the mass of the massive excitations can be computed in terms 
of correlators of the Heisenberg model. It is our purpose in this section to further explore this 
equivalence and to set up the formalism needed in what follows. 

In the thermodynamic limit the Schwinger Hamiltonian ( 2.33[ ), rescaled by the factor e\a/2, 
reads 

H = H Q + eH h (3.36) 

with 



N 



Ho = E 

x>y 

H h = -i(R - L) 



(x-yf 



(x - y) 



p{x)p{y) 



(3.37) 
(3.38) 



and e = t/e 2 L a 2 . In Eq.(3.38) the right R and left L hopping operators are defined (L = E)) as 



N 



N 2 



N 2 



x—l x—1 a—1 x— 

On a periodic chain the commutation relation 

[R, L]=0 

is satisfied. 



(3.39) 



(3.40) 
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We shall consider the strong coupling perturbative expansion where the Coulomb Hamiltonian 



(3.37) is the unperturbed Hamiltonian and the hopping Hamiltonian (3.38) the perturbation. Due 
to Eq. ( |2.34 ) every configuration with one particle per site has zero energy, so that the ground state 



of the Coulomb Hamiltonian ( |3.37 ) is 2 N times degenerate. The degeneracy of the ground state 



can be removed only at the second perturbative order since the first order is trivially zero. 

At the second order the lattice gauge theory is effectively described by the antiferromagnetic 
Heisenberg Hamiltonian. The vacuum energy — at order e 2 — reads 

4 2) =< Ht / H h > (3.41) 

where the expectation values are defined on the degenerate subspace and n is the operator projecting 
on a set orthogonal to the states with one particle per site. Due to the vanishing of the charge 
density on the ground states of Hq, the commutator 

[Ho, Hh] = Hh (3.42) 



holds on any linear combination of the degenerate ground states. Consequently, from Eq.( 3.41 ) one 
finds 

= -2 < RL > . (3.43) 

On the ground state the combination RL can be written in terms of the Heisenberg Hamiltonian. 
By introducing the Schwinger spin operators 



, _./,t ^± 

the Heisenberg Hamiltonian Hj reads 



S* = ^i, x -f^ (3.44) 



=i 

and, on the degenerate subspace, one has 

/ N 



x = l 
N 



\l x R x - ]p(x)p(x + 1) ) (3.45) 




<^>=(E ) = (E -5 iA ) • (3 - 46) 



Taking into account that products of L x and R v at different points have vanishing expectation 
values on the ground states, and using Eq.(|3~46|), Eq.(|3~43|) reads 



E { 2) = 4<Hj> . (3.47) 
The ground state of Hj singles out the correct vacuum, on which to perform the perturbative 



expansion. In one dimension Hj is exactly diagonalizable [|12j, 29 . In the spin model a flavor 1 
particle on a site can be represented by a spin up, a flavor 2 particle by a spin down. The spectrum of 
Hj exhibits 2 N eigenstates; among these, the spin singlet with lowest energy is the non degenerate 
ground state \g.s. >. 

We shall construct the strong coupling perturbation theory of the 2-flavor Schwinger model using 
\g.s. > as the unperturbed ground state. \g.s. > is invariant under translations by one lattice site, 
which amounts to invariance under discrete chiral transformations. As a consequence, at variance 
with the I-flavor model chiral symmetry cannot be spontaneously broken even in the infinite 
coupling limit. 

\g.s. > has zero charge density on each site and zero electric flux on each link 

p(x)\g.s.>=0 , E x \g.s. >= (x=l,...,N) . (3.48) 



7 



\g.s. > is a linear combination of all the possible states with ^ spins up and ^ spins down. The 
coefficients are not explicitly known for general N. In the appendix A we shall exhibit \g.s. > 
explicitly for finite size systems of 4, 6 and 8 sites. The Heisenberg energy of \g.s. > is known 
exactly and, in the thermodynamic limit, is [l4|] 

Hj\g.s. >= (-JV In 2)\g.s. > . (3.49) 



Eq.( |3.49|) provides the second order correction Eq.( 3.47 ) to the vacuum energy, Ef) s . = -42V In 2 



There exist two kinds of excitations created from \g.s. >; one kind involves only spin flipping 
and has lower energy since no electric flux is created, the other involves fermion transport besides 
spin flipping and thus has a higher energy. For the latter excitations the energy is proportional to 
the coupling times the length of the electric flux: the lowest energy is achieved when the fermion is 
transported by one lattice spacing. Of course only the excitations of the first kind can be mapped 
into states of the Heisenberg model. 

In p4| the antiferromagnetic Heisenberg model excitations have been classified. There it was 
shown that any excitation may be regarded as the scattering state of quasiparticles of spin-^: every 
physical state contains an even number of quasiparticles and the spectrum exhibits only integer spin 
states. The two simplest excitations of lowest energy in the thermodynamic limit are a triplet and 
a singlet fbif ; they have a dispersion relation depending on the momenta of the two quasiparticles. 
For vanishing total momentum (relative to the ground state momentum P g . s . = for ^ even, 
Pg.s. = 7T for ~ odd) in the thermodynamic limit they are degenerate with the ground state. 

In the appendix A we show that even for finite size systems, the excited states can be grouped 
in families corresponding to the classification given in pj. We explicitly exhibit all the energy 
eigenstates for N — 4 and N — 6. The lowest lying are a triplet and a singlet, respectively; they 
have a well defined relative (to the ground state) P-parity and G-parity — 1 h for the triplet and 
++ for the singlet. Since they share the same quantum numbers these states can be identified, in 
the limit of vanishing fermion mass, with the soliton-antisoliton excitations found by Coleman in 
his analysis of the 2-flavor Schwinger model. A related analysis about the parity of the lowest lying 
states in finite size Heisenberg chains, has been given in pp[ . 

Moreover in |l4| a whole class — M. af — of gapless excitations at zero momentum was singled out 
in the thermodynamic limit; these states are eigenstates of the total momentum and consequently 
have positive G-parity at zero momentum. The low-lying states of the Schwinger model also contain 
many massless excitations with positive G-parity; they are identified [^| with the excitations 
belonging to M.af- The mass of these states in the Schwinger model can be obtained from the 
differences between the excitation energies at zero momentum and the ground state energy. The 
energies of the states \ex. > belonging to the class Maf have the same perturbative expansion of 
the ground state. Consequently, the states \ex. > at zero momentum up to the second order in the 
strong coupling expansion have the same energy of the ground state ( 3.41 ), E\J = — 4iVln2. To 



this order the mass gap is zero. Higher order corrections may give a mass gap. 



4 The meson masses 

In this section we determine the masses for the states obtained by fermion transport of one site 
on the Heisenberg model ground st ate. Our analysis shows that besides the G-odd pseudoscalar 
isosinglet with mass mj = e-\/2/7r, there are also a G-even pseudoscalar isosinglet ++ and 
isotriplet 1 h and a G-odd scalar isotriplet 1 H with masses of the order of ms or greater. The 
quantum numbers are relative to those of the ground state I p f = ++ for N/2 even I PG = 0~~ 
for iV/2 odd. 

Two states can be created using the spatial component of the vector j 1 (x) Eq. ( p. 10 ) and isovector 



j\(x) Eqs.( ^8p^ ) Schwinger model currents. They are the G-odd pseudoscalar isosing let I PG 



and the G-even pseudoscalar isotriplet I PG = 1 + . The lattice operators with the correct 
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quantum numbers creating these states at zero momentum, when acting on \g.s. >, read 

N 

S = R + L = J2j\x) (4-50) 

x=l 

N 

(T_)t = +L< 12 ) = (4.51) 

(4.52) 



T+ = 



x=l 



V 2 



i?( afc ) and i> fc ) in ( gj^fjg ) are the right and left flavor changing hopping operators (L^ — 
(i?( afc ))t) 

JV 

R [ab) = T,<^ iA ^ ■ 

x=l 



The states are given by 



They are normalized as 

< S\S > 
<T+\T+ > 

and 



\S> = |0~ >= S\g.s. > 
|T±> = |1-+±1 >=T ± \g.s. > 
\T > = \l~+,0 >=T \g.s. > . 



< g.s.\SftS\g.s. >= -4 < g.s.\Hj\g.s. >= 4iVln2 
1(N+ < g.s.\Hj\g.s. >) = ^N(l - In 2) 



< T |T >=< T_|T_ >=< T + \T+ > 
In Eqs.(|4~5^,^57|,^5^) < g.s.\g.s. >= 1. 



(4.53) 
(4.54) 
(4.55) 



(4.56) 
(4.57) 

(4.58) 



The isosinglet energy, up to the second order in the strong coupling expansion, is E$ = E s 

„„ _ <S\Hq\S>_ 



(0) 



e 2 E { s 2) with 



< 5|5 > 



E 



,-, = < S\HlAsHh\S > 
< 5|5> 



As = (o? 5 — and 1 — lis a projection operator onto \S >. On Ip.s. > 

-Eg —-f?0 

[JT„, (ns^)"5] = (n + l)(n s F A ) n 5, (n = 0, 1, . 



holds; Eg. (4.60) may then be written in terms of spin correlators as 



4 2) = 4i+4- 



J2x=i < g- s -\ s x ■ s x+2 - \\g-s. > 
< g-s.\Hj\g.s. > 



(4.59) 
(4.60) 

(4.61) 
(4.62) 



One immediately recognizes that the excitation spectrum is determined once < g.s.\S x ■ S x +2\g-S- > 
is known. Equations similar to Eq.(4.62) may be established also at a generic order of the strong 
coupling expansion. 



E^ = 1. Following the same procedure as before one may compute the energy of the states ( 4.54 ) 



At the zeroth perturbative order the pseudoscalar triplet is degenerate with the isosing let E% 



(0) 
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and ( 4.55| ) to the second order in the strong coupling expansion. To this order, the energy is given 
by 



E, 



(2) 



Ei 2 } 



/CO 



4 < g.s.\Hj\g.s. > +5J2 x =i < gs.\S x ■ S x+2 - \\g.s. > 
N+ < g.s.\Hj\g.s. > 

where in terms of the vector operator V = ^2 x —i S x A S x+ i, one can write Ads(T) as 



A DS (T ± ) 
Ads (T ) 



= 12 
= 12 



< ff.s.|Q/i) 2 |.g.s. > + < g.s.\{V 2 ) 2 \g.s. > 

N+ < g.s.\Hj\g.s. > 
2<g. S .\(V 3 ) 2 \g.s.> 
N+ < g.s.\Hj\g.s. > ' 



(4.63) 

(4.64) 
(4.65) 



The v.e.v. of each squared component of V on the rotationally invariant singlet \g.s. > give the 
same contribution i.e. Aug(T±) ~ Aug(T ): the triplet states (as in the continuum theory) have a 
degenerate mass gap. This is easily verified by direct computation on finite size systems; when the 
size of the system is finite one may also show that A £>s is of zeroth order in N. 

2 2 

The excitation masses are given by ms — ^-{Es — E g . s .) and mj = ^y~(Et — E g , s ,). Conse- 

(2) (2) 

qucntly, the (iV-dependent) ground state energy terms appearing in E s and Ej, cancel and what 
is left are only N independent terms. This is a good check of our computation, being the mass an 
intensive quantity. 

In principle one should expect also excitations created acting on \g.s. > with the chiral currents, 
in analogy with the one flavor Schwinger model where, as shown in Ref. the chiral current 
creates a two-meson bound state. The chiral currents operators for the two flavor Schwinger model 
are given by 



The corresponding lattice operators at zero momentum are 



= *-£ = £x=ii 6 (*) 

= (T 5 ) t =jR ( 1 2)_ jL (12) =E W ii 5 (a;) 
UR(U) _ L (ll) _ R (22) + L (2 2)) = ^ if( 



T + 



V2 



The states created by (4.68 4.69,4.70) when acting on \g.s. >, are 

\S 5 > = |0++ >= S 5 \g.s. > 

|1+-,±1 >=Tl\g.s. > 
|T 5 > = |l+-,0>=T 5 | ff . S . > 



\Tl> = 



They are normalized as 

< S 5 \S 5 > 

< tI\tI > 

and 



< g.s.\S 5i <S 5 \g.s. >= -4 < g.s.\Hj\g.s. >= 47Vlog2 
^(N+ < g.s.\Hj\g.s. >) = ^N(l - log 2) 



< T°\T° >=< T_\T_ >=< Tl\Tl > 



Following the computational scheme used to study \S > and \T >, one finds for the state 



E 



E 



(o) 

(2) 
S 5 



^g.s. 



12-3 



Y,x=i < 9-s-\S x ■ S x+2 - \\g.s. > 
< 9-s-\Hj\g.s. > 



4.66) 
4.67) 

4.68) 
4.69) 
4.70) 
4.71) 

4.72) 
4.73) 
4.74) 

4.75) 
4.76) 

4.77) 

S 5 > 
4.78) 

4.79) 
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For the triplet |T 5 > one gets 

n(0) 



E. 



'I '5 



(2) _ 



1 



£(2) 
9-S- 



E I= i < ff-s-lffr ■ 5^+2 ~ |lg-£. > -4 < g.s.\Hj\g.s. > 
N+ < g.s.\Hj\g.s. > 



(4.80) 
(4.81) 



5 The correlators Ex < g.s.\S x ■ S x+ 2\g.s. > and < g.s.\V \g.s. > 

In this section we compute the spin-spin correlators needed to evaluate the second order energies of 
the isosinglet and isotriplets derived in section 4. The explicit computation of spin-spin correlation 
functions is far from being trivial since G(r) =< g.s.\So ■ S r \g.s. > is not known for arbitrary 
lattice separations r. For r = 2 it was computed by M. Takahashi in his perturbative analysis 
of the half filled Hubbard model in one dimension. For r > 2 no exact numerical values of G(r) 
are known. In jl^] were given two representations of G(r), while in J3l], |32| the exact asymptotic 
(r — > oo) expression of G{r) was derived. 



In order to explicitly compute the second order energies Eq.(4.62) and Eq.(4.63) one has to 
evaluate the correlation function 



1 - 

G(2) = -J2<9-s-\S x -S x+2 \g.s.> 

x=l 

which has been exactly computed in flj"7| and is given by 

G(2) = -(1 - 161n2 + 9C(3)) = 0.1820 



(5.82) 



(5.83) 



In the appendix B we shall show how the knowledge of this correlator allows one to compute 
explicitly the first three "emptiness formation probabilities", used in Ref.|18 in the study of the 
Heisenberg chain correlators, G{r). 

In order to compute the mass of the isotriplet one has to evaluate also the correlation functions 



appearing in Eq.(4.63); namely < g.s.\V -V\g.s. >, where V = J2 X S x AS x +i. The explicit expression 
of this function is not known. For its evaluation it is most useful to rearrange this correlator as 

N 

< g.s.\V ■ V\g.s. >= ]T (< g.s.\(S x ■ S y )(S x+1 ■ S y+1 )\g.s. > 

x,y=l 

N 

- < g.s.\(S x ■ S y+ i)(S x+ i ■ S y )\g.s. >) - J~] < g.s.\S x ■ S x+ i\g.s. > . (5.84) 



It is possible to extract a numerical value from Eq.(5.84) only within the random phase approxi- 
mation 0, ||. For this purpose it is first convenient to rewrite the unconstrained sum over the 
sites x and y as a sum where all the four spins involved in the v.e.v.'s lie on different sites, 

<g.s.\V-V\g.s.> = ]T (< g.s.\(S x ■ S V )(S X+1 ■ S 

y+l)\g- s - > 

y^x±l 

- 3 

- < g.s.\(S x ■ S y+1 ){S x+1 ■ S y )\g.s. >) + -N 

^ N N 

- - 22 < g-s.\S x ■ S x+ i\g.s. > ~22 < 9-s.\S x ■ S x+2 \g.s. > (5.85) 



and then factorizc the four spin operators in Eq.(5.85) as 
< g.s.\V -V\g.s. > 



N^(< g.s.\S ■ S r \g.s. > 2 - < g.s.\S ■ S r+ i\g.s. >< g.s.\Si ■ S r \g.s. >) 

r=1 
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N N 



3 1 J - - v— - -. 

a ^ ~~ 9 X! < 9- s -\ S x ■ S x+1 \g.s. > - ^ < g.s.^ ■ 5 x+2 |3.s. > .(5.86) 



=i 



Of course, Eq.( |5.86 ) provides an answer larger than the exact result; terms such as < (. . .)(. . .) > 



contain negative contributions which are eliminated once one factorizes them in the form < (...) >< 
(...) >. This is easily checked also by direct computation on finite size systems. 

The spin-spin correlation functions G(r) are exactly known for r = 1,2. For the spin-spin 
correlation functions G{r) up to a distance of r = 30 the results are reported in table (pi) |3l], p4|. 



Table 1: Spin-spin correlation functions 



r 


G(r) 


r 


G(r) 


1 


-0.4431 


16 


0.0305 


2 


0.1821 


17 


-0.0296 


3 


-0.1510 


18 


0.0274 


4 


0.1038 


19 


-0.0267 


5 


-0.0925 


20 


0.0249 


6 


0.0731 


21 


-0.0242 


7 


-0.0671 


22 


0.0228 


8 


0.0567 


23 


-0.0223 


9 


-0.0532 


24 


0.0211 


10 


0.0465 


25 


-0.0206 


11 


-0.0442 


26 


0.0196 


12 


0.0395 


27 


-0.0193 


13 


-0.0379 


28 


0.0183 


14 


0.0344 


29 


-0.0181 


15 


-0.0332 


30 


0.0172 



For r > 30, one may write |3] 



G(r) 



with g(r) satisfying 



and 



3/2 1 



[1 - ^9{rf + 156C S ? V ) 3 + 0(g(rf) 



4 V 7T 3 r^fgjr) 16 384 

^((-l)r + l + O(0(r)) + O(i)] 
2r H 



9(r) = -^T-r{l + ^g(r)ln(g(r))) 



C{r) = ln(2V2^e 7+ V) . 
Eq.(5.88) may be solved by iteration. To the lowest order in ^ one finds 



Inserting (5.90) in (5.87) leads to 



G(r) « V2^VC(r)[l + ^ lnC(r)] + O(^) 



(5.87) 

(5.88) 
(5.89) 

(5.90) 
(5.91) 
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Putting ( pT9lj ) in ( pT86| ), one finally gets 

< ff.s.|y.?|s.s. >=0.3816iV 



Using Eq.( |5.83[ ), the isosinglet mass reads as 



7719 1 t 

-5 - = - + 1.9509 e 2 

e z a 2 



For what concerns the isotriplet mass, since the double sum in Eq.( 4.63] ) is given by 

< g.s.\V -V\g.s. > 



Ads(T) = 8 



N+ < g.s.\Hj\g.s. > 



using Eq. fl5.92D , one gets 



rriT 

2 

ei a 



- + 0.0972 e 2 



(5.92) 



(5.93) 



(5.94) 



(5.95) 



The existence of massive isotriplets was already noticed in pi| , and their mass in the continuum 
theory was numerically computed for various values of the fermion mass. In particular there is a 
G-parity even isotriplet with mass approximately equal to the mass of the isosinglet . 

The mass of the I1S5 > isosinglet and the \T$ > isotriplet is 



7775.5 



777 T 5 



1 9 

= - + 5.85e 2 
= - + 4.4069e 2 



(5.96) 
(5.97) 



Equations ( 5.93| ), (5.95), ( 5.96| and (5.97) provide the values of ms, rnx and 777,7-5 for small 
values of z — e 2 — U P to the second order in the strong coupling expansion. Whereas ( 5.95| ) is 

only approximate (5.93) and ( |5.97 ) are exact at the second order in the e expansion. In section 6 
we shall extrapolate these masses to the continuum limit using the standard technique of the Pade 
approximants. 



6 The chiral condensate 

In the following we shall first prove that also on the lattice either the isoscalar Otpip) or the isovector 
(tjj<j a ip') chiral condensates are zero to every order of perturbation theory. This should be verified 
by explicit computation since on the lattice the symmetry SUl(2) ® SUr(2) is already broken by 
introducing staggered fermions; thus, there is no symmetry to prevent the formation of such chiral 
condensates. In the continuum theory, instead, the breaking of the SUl(2) ® SUr(2) down to 
SUv{2) is prevented by the Coleman theorem Jl9| . 

In the staggered fermion formalism the isoscalar condensate is given by 

Y / Mx)M^^ { -^(4,^i, x + 4,^2, x -i) ; (6.98) 

a— 1 

it is obtained by considering the mass operator 

1 N 

M = 2N~ a Et-^W'LV'i,* + 4,M (6.99) 

x-l 

and evaluating its expectation value on the perturbed states \p g . s . > generated by applying Hh to 
\g.s. >. To the second order in the strong coupling expansion, one has 

Ws. >= \g-s. > +e\pl s . > +e 2 \ P 2 g . s . >+... (6.100) 
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where 

\p l g .s. > = -H h \g.s. > (6.101) 
\p 2 g ,. > = ^H h H h \g.s. ' ■ (6-102) 



To the fourth order, ( 6.99 ) is given by 

_ <P 9 .s.\M\ Pg . s . > _ < as. M n.s. > +e 2 < pl. s .\M\ P l s . > +e 4 < p 2 g . s \M\p 2 g . s . > 



X. 



< Pg. s .\Pg.s. > < g-s.\g-s. > +e 2 < pl. s .\pl. s . > +e 4 < p 2 . s .\p 2 . s . > + ... 

(6.103) 

It is very easy to see that Xisos. is zero to all orders in the strong coupling expansion. Let us 
introduce the translation operator 

T = e lj3a ; (6.104) 



using 



and 



TMf- 1 = -M (6.105) 
THhf- 1 = H h (6.106) 



one gets order by order in the strong coupling expansion in Eq.(6.103) 



T\g.s.>=±\g.s.> (6.107) 

(6.108) 



In Eq.(6.107) the + appears when N/2 is even and the — when N/2 is odd. 

The isovector chiral condensate is given by the expectation value of the operator 

i N 

^ = ^N~ a Ec-wU^.* ( 6 - 109 ) 

x=l 

on the perturbed states \p g . s . >■ Taking into account that 

TZf- 1 = -S (6.110) 

one gets 

/Cisov. — Xisov. i (6.111) 

also the isovector chiral condensate is identically zero. 

In the continuum there is, as evidenced in section 2, only a non-vanishing chiral condensate 
associated to the anomalous breaking of the Ua(^) symmetry |24|, |2q| . Since 

^ L ( x )^ R ( x )=^( x )l±pLr{x) , (6.H2) 

the pertinent operator is F = 4^ ^r^ 4^r \ its expectation value has been computed in pH ]25|] 
and is given by 

< ^ViMW >= (£) 2 ^ = (£) 2 m| . (6.113) 
On the lattice one has 

^ i (x)^( SB )^^i(^ 0lX -^, B+1 V'-,^i+4 a) -4 0) ) • (6-114) 

The factor l/2a is due to the doubling of the lattice spacing in the antiferromagnetic bipartite 
lattice. Upon introducing the occupation number operators = x ip a ,xi t ne umklapp operator 
F is represented on the lattice by 

F = ^ N { (n ^ ~ n *+ l)(n * ) " n ™J + {L * } ~ 41))(42) " ^ 2)) } • (6 - 115) 

X — l 
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The strong coupling expansion carried up to the second order in e = , yields 



<Fy _ <p g . s .\F\p g . s . > ■ .'/•-• //•■-• • i>;,.. > (6n6) 

<Pg.s.\Pg.s.> < g-S.\g.S. > +£ 2 < p 1 ]p X , > 



for the lattice chiral condensate. Since 



<g.s.\g.s.> = 1 (6.117) 

< Pistils. > = -A<g.s.\Hj\g.s.> , (6.118) 

and taking into account that 

<g.s.\F\g.s.> = < g.s.\Hj\g.s. > (6.119) 
<pIs.\F\pIs.> = -^(-2<g.s.\(Hj) 2 \g.s.>~<g.s.\Hj\g. S .>+^N 

2 - 1 

~ 3E <g-s.\S x -S x+2 --\g.s, >) , (6.120) 

X—l 



from Eqs.(|3.49|) and ( |8.223| ), one gets 

< F >= 4-(0-0866-0.4043e 2 ) . (6.121) 

or 

A nonvanishing value of the lattice chiral condensate is due to the coupling — induced by the lattice 
gauge held — between the right and left fermions. This is the relic in the lattice of the Ua{1) 
anomaly in the continuum theory. 

7 Lattice versus continuum 

We now want to compare our lattice results with the exact results of the continuum model; to do 
this, one should extrapolate the strong-coupling expansion derived under the assumption that the 
parameter z — e = 4 4 <C 1 to the region in which z 3> 1; this corresponds to take the continuum 

limit since e L a A — ► when z — > 00. To make the extrapolation possible, it is customary to make 
use of Pade approximants, which allow to extrapolate a series expansion beyond the convergence 
radius. Strong-coupling perturbation theory improved by Pade approximants should then provide 
results consistent with the continuum theory. As we shall see the strong-coupling expansion derived 
in this paper provides accurate estimates of the meson masses, already at the hrst order in powers 
of z. 

Let us now evaluate ms and the lattice light velocity t. We hrst compute the ratio between the 
continuum value of the meson mass ms = \J^ e c and the lattice coupling constant by equating 
the lattice chiral condensate, Eq.( |6.121 ), to its continuum counterpart Eq.( 6.113| ) 



1 ^ 

^■(0.0866 - 0.4043z) = (^—fml . (7.122) 



Eq.(7.122) is true only when Pade approximants are used since, as it stands, the left hand side 
holds only for z <C 1, while the right-hand side provides the value of the chiral condensate to be 
obtained when z — 00. Using 

1 

a=-^ , (7.123) 



one gets from Eq.( 7.122| ) 



, ms \ 1 , 47r „ , z 2 , 

— ) 2 = — ) 2 — 0.0866 - 0.4043z . (7.124 
eh e 7 t 
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As in Refs.|3|, [|, due to the factor za, the second power of Eq.(7.124) should be considered in 
order to construct a non diagonal Pade approximant. Since the strong coupling expansion has been 
carried out up to second order in z, one is allowed to construct only the [0, 1] Pade approximant 
for the polynomial written in Eq.(7.124). One gets 



,m s )4 _ ( 4tt 4 1 0.0074Z 
[ e L ' ~ W t 2 1 + 9.3371z 



and, taking the continuum limit z — > oo, one finds 



m s 4 4tt 4 0.0008 



e L V 



t 2 



(7.125) 



(7.126) 



Next we compute the same mass ratio by equating the singlet mass gap given in Eq.(5.93) to 
its continuum counterpart ms 



9 A 

e\a{- + 1.9509z) 



m s 



(7.127) 



Again, Eq.(7.127) is true only when Pade approximants are used. Dividing both sides of Eq.( |7.127j ) 
by e_L and taking into account that 



one gets 



e L a = — 

Zi 



22 = 4.(1 + 1.9509*) 

eL zi 2 



(7.128) 



(7.129) 



Taking the fourth power and constructing the [1,0] Pade approximant for the right hand side of 
Eq.( |7.129| ) one has 

t 2 , 1 



,ms_ A 
eL 



-(- + 0.9754z) ; 
z 16 



when z — > oo, Eq.(7.130) gives 



»4 =t 2 9754 

eL 



(7.130) 
(7.131) 



The numerical value of the hopping parameter t, determined if one equates Eq.(7.131) and 
Eq.( |7.126| ), is 

4-7T 

t = —0.1692 = 1.1940 (7.132) 



and lies 19% above the exact value. Putting this value of t in Eq.( 7.126j ) or Eq.(7.131) one gets 

m s 



= 1.0969 



PL 



(7.133) 



which lies 37% above the exact value y ^- It is comforting to see that the lattice reproduces in a 
sensible way the continuum results even if we use just first order (in z) results of the strong coupling 
perturbation theory. 

Using the value of t given in Eq.( [7.132 ) one gets for the isotriplet mass Eq.(7.139) 



eL 



0.5143 



By direct computation on an 8 sites chain one gets 

eL 



1.3524 



(7.134) 



(7.135) 



The discrepancy between Eq.(7.134) and Eq.(7.135) is mainly due to the approximation involved 
in the computation of < g.s.\V 2 \g.s. >. However, it is safe to believe that our lattice computation 
implies the existence of a massive isotriplet 1 v with a mass between the lower bound (7.134) and 
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the upper bound ( 7.135| ). This is in agreement with the results provided for the continuum theory 

in gg. 

Using a similar procedure one may also compute the masses of the singlet ++ and the triplet 
From Eqs.( |5.96||T97| ) one gets 

— = 1.4290 . (7.136) 
= 1.3347 . (7.137) 



The triplet \T 5 > , being G-odd, is a scattering state of a singlet with a 1 h triplet, which are 
the fundamental excitations of the system. The mass of this 1 H triplet should be larger than the 
mass of the massive 1 v triplet, which should be a scattering state of massless 1 h triplets. 

Putting t = c = 1, i.e. e^a = -\- the lattice mass spectrum gets closer to its continuum 

counterpart; for the isosinglet mass, one gets 

— = 0.9938 (7.138) 

e 

while for the isotriplet one gets 

— = 0.4695 . (7.139) 

eL 

Eq.( 7.138 ) provides a value of the isosinglet mass lying 24% above the exact answer. Again the 
triplet mass is reproduced with lesser accuracy due to the random phase approximation used in the 
computation of the pertinent correlator; a better answer is given however by a direct computation 
on the 8 sites chain yielding the value 1.2346 for rriT/eL- 



8 Concluding remarks 

In this paper we used the correspondence between the 2-fLavor strongly coupled lattice Schwinger 
model and the antiferromagnetic Heisenberg Hamiltonian established in Q to investigate the spec- 
trum of the gauge model. Using the analysis of the excitations of the finite size chains given in the 
appendix, we showed the equality of the quantum numbers of the states of the Heisenberg model 
and the low-lying excitations of the 2-flavor Schwinger model. We provided also the spectrum of 
the massive excitations of the gauge model; in order to extract numerical values for the masses, 
we explicitly computed the pertinent spin-spin correlators of the Heisenberg chain. Although the 
spectrum is determined only up to the second order in the strong coupling expansion the agreement 
with the continuum theory is satisfactory. 

The massless and the massive excitations of the gauge model are created from the spin chain 
ground state with two very different mechanisms: massless excitations involve only spin flipping 
while massive excitations are created by fermion transport besides spin flipping and do not belong to 
the spin chain spectrum. As in the continuum theory, due to the Coleman theorem f20|j , the massless 
excitations are not Goldstone bosons, but may be regarded as the gapless quantum excitations of 
the spin-i antiferromagnetic Heisenberg chain fl35| . 

In computing the chiral condensate we showed that, also in the lattice theory, the expectation 
value of the umklapp operator F is different from zero, while both < tpij) > and < ipa a ?p > are 
zero to every order in the strong coupling expansion. This implies that both on the lattice and the 
continuum the SU(2) flavor symmetry is preserved whereas the Ua(1) axial symmetry is broken. 
The umklapp operator F is the order parameter for this symmetry, but being quadri-linear in the 
fermi fields, is invariant, in the continuum, under chiral rotation of tt/2 and on the lattice under the 
corresponding discrete axial symmetry ( 2.30D (translation by one lattice site). This shows that the 



discrete axial symmetry is not broken in both cases. Our lattice computation enhance this result 
since the ground state of the strongly coupled 2-flavor Schwinger model is translationally invariant. 

The pattern of symmetry breaking of the continuum is exactly reproduced even if the Coleman 
theorem does not apply on the lattice and the anomalous symmetry breaking is impossible due to 
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the Nielsen-Ninomiya [[38| theorem. At variance with the strongly coupled 1-flavor lattice Schwinger 
model, the anomaly is not realized in the lattice theory via the spontaneous breaking of a residual 
chiral symmetry [Q, but, rather, by explicit breaking of the chiral symmetry due to staggered 
fcrmions. The non-vanishing of < F > may be regarded as the only relic, in the strongly coupled 
lattice theory, of the anomaly of the continuum 2-fiavor Schwinger model. It is due to the coupling 
induced by the gauge field, between the right and left-movers on the lattice. 

It is quite straightforward to generalize our analysis to an SU {A/")-flavor group. For this, one 
should observe || that the results are very different depending on if Af is even or odd. When Af is 
odd, the ground state energy — in the strong coupling limit — is proportional to e 2 but, when Af is 
even, the ground state energy is of order 1. This difference arises since, on the lattice, the charge 
density operator is odd under charge conjugation; therefore, the constant 4r should be subtracted 
from the charge density operator 0. As a consequence, when Af is odd, the ground state supports 
electric fluxes while this becomes impossible when Af is even. 

When the fermion mass m is different from zero, some further difference arises between Af odd 
and Af even. When Af is odd, the mass term induces a translational non invariant ground state, 
generating a spontaneous chiral symmetry breaking. When Af is even, the ground state remains 
translationally invariant in the strong coupling limit, i.e. e 2 S> m. In the weak coupling limit, 
m ^> e 2 , the discrete chiral symmetry is spontaneously broken for every Af. For Af = 2, the 
soliton-antisoliton excitations f2C(| acquire a mass. 
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Appendix A: 

Exact diagonalization of finite size antiferromagnetic Heisen- 
berg chains 

In the following we shall provide the exact diagonalization of the Heisenberg antiferromagnetic 
model for finite size chains of N — 4, 6 and 8 sites and compare the results with the Bethe Ansatz 
solution provided in fll4| . In order to make our arguments self-contained, we shall outline the steps 
involved in deriving the exact solution in the thermodynamic limit. We shall show that already very 
small finite size chains exhibit spectra that match very well with the thermodynamic limit solution. 
Furthermore the analysis of finite size chains is very useful since it allows the comparison between 
the quantum numbers of the Schwinger and the Heisenberg model excitations. The interested reader 
may also look up references |l4|, |3(| |29| . 

The one dimensional isotropic Heisenberg model describes a system of N interacting spin-i 
particles. The Hamiltonian of the model is 

N 1 
Hj = J^2(S x -S x+l - -) . (8.140) 

x—l 

where J > ( J < would describe a ferromagnet) and the spin operators have the following form 

S x = li <g> 1 2 <8 y ® ...<8> ljv ■ (8.141) 
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They act nontrivially only on the Hilbert space of the X th site. Periodic boundary conditions are 
assumed. 



The Hamiltonian Q3.140 ) is invariant under global rotations in the spin space, generated by 

JV 

S = ^§ x . (8.142) 

x = l 

Due to the periodic boundary conditions, under translations generated by the operator T, 

f : S x S x -! (8.143) 

the Hamiltonian is invariant and [S, T] = 0. 

The energy and the momentum of a given state with M spins down can be written as fl2]| 

M M 

E M = ^ £q = -_^ (8.144) 

a—l a—1 a 4 



M M , 



A I 



ilnT=Y,p a =iJ2^^- r ■ (8-145) 



Act ~t~ o 
a=l a=l 



Energy and momentum are thus additive as if there were M independent particles and the param- 
eters A Q must satisfy the Bethe Ansatz equations 

M^-IIr^fe . ( 8 - 146 ) 

in order for and P m to b e eigenvalues of the Hamiltonian and momentum operator. See fl4f| 
for a derivation of Eqs.( ^.14^ ). 

The solution of the antiferromagnetic Heisenberg chain is reduced to the solution of the system 



of the M algebraic equations ( S.146 ). This, in general, is not an easy task. It can be shown 14 



however, that, in the thermodynamic limit N —> oo, the complex parameters A have the form 

Xl = Xj,L+H , 1= -L,-L + 1,...,L- 1,L; (8.147) 



where L is a non-negative integer or half- integer, A^l is the real part of the solution of ( 8.146| ) 
and we shall define shortly the set of allowed values for the integer index j. The A's that, for a 
given Ajr,, are obtained varying I between [-L, L] by integer steps, form a string of length 2L + 1, 
see fig.(p. This arrangement of A's in the complex plane is called the "string hypothesis" [[u]. In 
the following we shall verify that, even on finite size systems, the "string hypothesis" is very well 
fulfilled. 



In a generic Bethe state with M spins down, there are M solutions to ( S.146) ), which can be 



grouped according to the length of their strings. Let us denote by vl the number of strings of 
length 2L + 1, L = , i, . . .; strings of the same length are obtained by changing the real parts, Aj.L, 



of the A's in ( p.l47| ); as a consequence j = 1, . . . , vi,. If one denotes the total number of strings 
by q one has 

q = Y. v L > M = J2( 2L + 1 >L ■ (8.148) 



The set of integers (M, q, {vl}) constrained by ( 8.148| ), characterizes Bethe states up to the 



fixing of the q numbers \j l\ this set is called the "configuration". Varying M, q and v^, one is 
able to construct all the 2 N eigenstates of an Heisenberg antiferromagnetic chain of N sites |14 



The energy and momentum of the Bethe's state, corresponding to a given configuration — within 
exponential accuracy as N — > oo — consist of q summands representing the energy and momentum 
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Figure 1: Strings for L = 0, | 



of separate strings. For the parameters Xj.l of the given configuration, taking the logarithm of 
(8.146) the following system of equations is obtained in the thermodynamic limit 



2N arctan ■ 



Li 



hi k=l 



where 



®l 1 l 2 W = 2 ^ (arctan ^ + arc tan 
L=|L 1 -L 2 |^0 



£ + 1' 



(8.149) 



.150) 



Integer and half integer numbers Q^l parametrize the branches of the arcotangents and, conse- 
quently, the possible solutions of the system of Eqs.( ^.149 ). In ref.[jL4j it was shown that the Qj,t 
are limited as 

- < Q hL < Q 2 + <...< Q Vl ,l < 



ymax 



with Q™ aa: given by 



and 



yinax 



N 
~2 



YJ(l,l')v l ,-1 



J(L U L 2 



2min(Li, L 2 ) + 1 if L 1 ^ L 2 
2L t + \ if Lx = L 2 



(8.151) 
(8.152) 

(8.153) 



The admissible values for the numbers Qj,L are called the "vacancies" and their number for every 
L is denoted by Pl 

Pl = 2Q" L iax + 1 . (8.154) 

The main hypothesis formulated in (m) is that to every admissible collection of Qj,t there 
corresponds a unique solution of the system of equations ( 3. 149 ) . The solution always provides, in 
a multiplet, the state with the highest value of the third spin component S 3 . 

Let us now consider some simple example. The simplest configuration has only strings of length 
1, i.e. all the A's are real. The singlet associated to this configuration 

N 

M = q = v = — , v L = , L>0 , (8.155) 
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is the ground state. The vacancies of the strings of length 1, i.e. the admissible values of Qj t o, due 
to eqs.( fe,151 8.152,8.153), belong to the segment 



N 1 N 1 

h-<Qin< 

4 2- - 4 2 



.156) 



Therefore they are A/2. All these vacancies must then be used to find the A/2 strings of length 1. 
As a consequence this state is uniquely specified and no degeneracy is possible. 

Next we consider the configuration that provides a singlet with 1 string of length 2 and all the 
others of length 1: 



2 



2 



Vn = 2 

2 



1 



= 



r 1 

L> 2 



.157) 



For the strings of length 1 the number of vacancies is again N/2; for the string of length 2 there is 
one vacancy and the only admissible Qj t i equals 0. Thus, since the number of strings of length 1 is 
^0 = ^ — 2, there are two vacancies for which Eqs.( 8.149| ) have no solution; they are called "holes" 



and are denoted Qi and ■ This configuration is det ermin ed by two parameters: the positions 



of two "holes" which vary independently in the interval ( 8.156 ). 

There is another state with only 2 holes: the triplet corresponding to the configuration 



M 



N 
~2 



1 



L > 



(8.158) 



The number of vacancies for the strings of length 1 equals + 1 , while — ^ — 1 



The excitations determined by the configurations ( 8.157 , 8.158 ) belong to the configuration class 
called in jlj] Mat- The class Mat is characterized as follows: the number of strings of length 
1 in each configuration belonging to this class differs by a finite quantity from A/2, v§ — — fco 
where fco is a positive finite constant, so that the number of strings of length greater than 1 is finite. 
From (8.154) we then have 



A 



= - + fc -2£ 



VL 



L>0 



P L = 2k -2^2 J(L, L')v v , L > 



L'>0 



so that 



Po> 



A 



P L < 2k , L > 



(8.159) 
(8.160) 

(8.161) 



From (8.159) follows that the number of holes for the strings of length 1 is always even and equals 



2 only for the singlet and the triplet excitations discussed above. One can imagine the class M.af 
as a "sea" of strings of length 1 with a finite number of strings of length greater than 1 immersed 
into it. It was proven in |l4| that, in the thermodynamic limit, the states belonging to M. at have 
finite energy and momentum with respect to the antiferromagnetic vacuum, whereas each of the 
states which corresponds to a configuration not included in the class M. at has an infinite energy 
relative to the antiferromagnetic vacuum. 



Let us now sketch the computation of the thermodynamic limit ground state energy. Eqs.(B.14£) 
for the ground state have the form 



ttQ- 1 N/2 
arctan2Aj = ^ J + — arctan(Aj — Afc) 



fc=i 



Taking the thermodynamic limit N — > oo, one has 



N 



1 1 

7 < ^ < 7 
4 4 



(8.162) 



(8.163) 
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and Eqs.(8.162) can be rewritten in the form 

arctan2A(x) = irx + I arctan(A(a;) — \{y))dy 



Upon introducing the density of the numbers X(x) in the interval dX 



d\(x) 



*{X) I 

dx W=x{\) 



and differentiating Eqs.(8.164), one gets 



P(A) 



e 2 



2vr / l + e -|«l 



-iA|S| 



i 



2 cosh7rA 



(8.164) 



(8.165) 



(8.166) 



The density p{\) introduced in this way is normalized to 1/2. It is now easy to compute the energy 
and the momentum of the ground state 



2 /.OO 

Eg. s . =J2^ = N e(X)p(X)dX 

N_ 

P g .s. = Y.P- = N P(A)p(A)dA 

„ 1 J — OO 



JN 
"~4~ 

4/ dA- 



1 



(A 2 + \) cosh7rA 



N 



= -JN\n2 



cosh 7rA 2 



— 7r (mod 2tt) 



.167) 



(8.168) 



According to Eq. fl8.168D , P g . s . = (mod27r) for f even, and P 5 . s . = tt (mod 2?r) for f odd. The 
ground state, as expected, is a singlet, in fact the spin S is given by 



N ^ N 



S = — - ^l = ll-iV / j o(A)dA = 



2 ^ 2 



(8.169) 



Let us analyze the triplet described by Eq. fl3.158 ); Eqs.(8.149) take the form 



arctan2A 7 = ^ H arctan(A, — A^) 

J AT M Z-^i ■> 



N N 



.170) 



k=l 



where now the numbers Qj lie in the segment [— -j, ^] and have two holes, Q± and Q^ 1 ' with 
Qi < Q2 ■ Taking the thermodynamic limit one gets 



X 2 



N 1 N 

where 9{x) is the Heaviside function. Eqs.( 8~l~70|) become 



arctan2A(ir) = -ex + — {6{x — x\) + 8(x — X2)) + / arctan(A(x) — X(y))dy 



Eq. (|S.172| ) gives, for this triplet, the density of A, p(\) = ^ 



Pt (\) = p(X) + -HA - Ai) - a(X - A 2 )) 



where p(A) is given in Eq.( ^.166| ) and 



Z7T 



1 



l + e-l«l 



e^d^ 



.171) 



.172) 



.173) 



.174) 
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Ai and A2 are the parameters of the holes, A,; = A(#i), i = 1,2. The energy and the momentum of 
this state measured from the ground state are now easily computed 



(8.175) 
(8.176) 



e r (Ai, X 2 )=N e(A)(p t (A) - p(X))d\ = e(A x ) + e(A 2 ) 

J —OO 

Pr(Ai,A 2 ) =N / p(A)(pt(A)-p(A))dA = p(Ai)+p(A a ) (mod 2tt) 



where 



e(A) 



e(p)a(X — p)<i/i = J 



2 cosh 7rA 



P(A) 



p(fj,)a(\ ~ [ijdfi — arctan sinh 7rA - 



-7T < p(A) < 



From Eqs.( |8.177| , p.l78D one gets 



Jtt 



■ smp 



(8.177) 
(8.178) 

(8.179) 



The momentum pj^Ai, A2) varies over the interval [0, 2ir), when Ai and A2 run independently over 
the whole real axis. The spin of this state can be computed by the formula 



S = - 



(<r(A- A x ) + a(A- X 2 ))dX = 1 



(8.180) 



Let us finally analize the singlet excitation characterized by the configuration ( 3.157). Denoting 
by As the only number among the A^i/2 which characterizes the string of length 2 and by Xj the 
numbers A^o for the strings of length 1, Eqs.(3.149) read 



arctan 2Aj 
arctan A 5 



7TC 



1 



N N 
N 



^(Xj - A s ) + — ^2 arctan(Aj - X k 



fe=i 



I E - Xj) 



.181) 



(8.182) 



with 



$(A) — arctan 2 A + arctan - A 

3 



The — 2 numbers Qj vary in the segment [— 4 ., . , 



(8.183) 

o 

1 N — i]; among them there are the two holes 



'3 ""J """" u.^^.v L 4 ' 2' 4 2 

and Q^ 1 - Taking the thermodynamic limit one finds the density of A's for the singlet 



p(X s ) = p(X) + ^(ff(A - Ai) + a(X - A 2 ) + w(A - A s )) 



where p and cr were given in Eqs.( |B.166| , |8.174| ) and where 



w(A) 



1 

2^ 



-i|f|-iA« 



e 2 



tt(1 + 4A 2 ) 



(8.184) 



(8.185) 



In 1 14 it was demonstrated that the string parameter As is uniquely determined by the A's 
parametrizing the two holes 



As = 



A^ + A^ 



(8.186) 



In 113] it was also proved the remarkable fact that the string of length 2 does not contribute to the 
energy and momentum of the excitation, so that the singlet and the triplet have the same dispersion 
relations 



es(Ai,A 2 ) = £t(Ai, A 2 ) = e(Ai) + e(A 2 ) 

Ps(Ai,A 2 ) = Pt(Ai, A 2 ) = p(Ai) +p(A 2 ) (mod 2ir) 



(8.187) 
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The spin of this excitation is, of course, zero 



S 



(2cr(A) +uj(\))d\ = 



(8.189) 



.158) and the state of 



The only difference between the state whose configuration is given in Eq. 
Eq. flS. 1571 ) is the spin. 

To summarize, the finite energy excitations of the antiferromagnetic Heisenberg chain are only 
those belonging to the class M.at and are described by scattering states of an even number of 
quasiparticlcs or kinks. The momentum p of these kinks runs over half the Brillouin zone — tt < 
p < 0, the dispersion relation is e(p) = 4p sinp, Eq.(B.179), and the spin of a kink is 1/2. The 
singlet and the triplet excitations described above are the only states composed of two kinks, the 
spins of the kinks being parallel for the triplet and antiparallel for the singlet. For vanishing total 
momentum all the states belonging to Mat have the same energy of the ground state so that they 
are gapless excitations. Since the eigenstates of Hj always contain an even number of kinks, the 
dispersion relation is determined by a a set of two-parameters: the momenta of the even number of 
kinks whose scattering provides the excitation. There are no bound states of kinks. 

Let us now turn to the computation of the spectrum of finite size quantum antiferromagnetic 
chains by exact diagonalization. We shall see that already for very small chains, the spectrum is 
well described by the Bethe ansatz solution in the thermodynamic limit. Furthermore, an intuitive 
picture of the ground state and of the lowest lying excitations of the strongly coupled 2-flavor lattice 
Schwinger model emerges. 

The states of an antiferromagnetic chain are classified according to the quantum numbers of 
spin, third spin component, energy and momentum \S, S s ,E,p >. For a 4 site chain the momenta 
allowed for the states are: 0, ^ mod 2ir. The ground state is 



\g.s. >= |0, 0, -3J, >= 



1 



(2| UU> +2| 4T4T> -| TTU> -I TUT> -| UH> -| 4TTI>) 



(8.190) 
the P-parity 
> appearing 



This state is P-parit y even. In fact, by the definition of P-parity given in Eq.( 2.2£ ) 
inverted state ( 8.190 ) is obtained by reverting the order of the spins in each vector | . . 

in (jugg ), e.g. | UTT>^ | TTU>. 

The A's associated to the ground state (solution of the Bethe ansatz equations (3.146)) are 
Ai = K= and A 2 = —^=. There is also an excited singlet 

2v3 2V3 



|0,0,-J,7T >= 



7t il UTT> 



-UTU>-|TUT>+I TTU>) 



.191) 



It is P-even, so that it is a S p = + excitation, with the same quantum numbers (the isospin is 
replaced by the spin) of the lowest lying s inglet excitation of the strongly coupled Schwinger model 
discussed by Coleman 
by the configuration 



0|. The state (8.191) also coincides with the excited singlet described 



.157). 



approximately of length 2, A] 



It has only two complex A's which arrange themselves in a string 
48 g~ 17 and there are two holes with = 



-A, 



and 



There are also three excited triplets, whose highest weight states are 



IM, 

11,1,-2./.;: 
IM, 



7 37F 

J 'T > 



-4(1 ittt> +i\ Utt> -I TUt> -»| TTU>) 

V4 

-4(1 ittt> -I TITT> +1 TUT> -I TTTI>) 

V4 

4(1 ittt> -i\ TITT> -I TUT> +»| TTU>) 

V4 



(8.192) 
(8.193) 
(8.194) 



Among these, only the non-degenerate state with the lowest energy has a well defined P-parity 
(8.193). It is a S p = 1~ like the lowest lying triplet of the 2-flavor strongly coupled Schwinger 
model. The degenerate states can be always combined to form a P-odd state. 
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We thus see that within the states in a given configuration there is always a representative state 
with well defined parity, the others are degenerate and can be used to construct state of well defined 
energy and parity. Moreover the parity of the representative states (with respect to the parity of 
the ground state) is the same of the one of the lowest-lying Schwinger model excitations in strong 
coupling. 



All the triplets in (3.194) have one real A and two holes; they can be associated with the family 
of triplets ( p. 158 ). In table (|J) we summarize the triplet A's and Q^'s. The spectrum exhibits also 



Table 2: Triplet internal quantum numbers 



TRIPLET 


A 




|l,l,-J,f > 


l 

2 


-1 


|1,1,-2J,7T> 





-l i 


|M,-J,tF > 


1 

2 


1 



a quintet, whose highest weight state is 

|2,2,o,o >= | TTTT> 

We report in fig.(@) the spectrum of the 4 sites chain. 



E/J 



[3- 



tt/2 



3tt/2 2tt 



& B ^ 

© 



■a 



O QUINTET 
Q TRIPLETS 
□ SINGLETS 



P 



(8.195) 



Figure 2: Four sites chain spectrum 

Let us analize the spectrum of the 6 site antiferromagnetic chain. The momenta allowed for the 
states are now 0, |, ^r, n, ^ mod 27r. The ground state is 

!<?•*•> = |0,0,-^(5 + Vl3),7r>=-=i={|mUT>-|UTm> 
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l - VT3 

6 



(I mm> -i umt> +i mm> -i mm> +1 umt> -i mm> 



- I I4T4TT> +1 mtU> -I UTTU> +1 iTTUT> -| TTiiti> +1 TUUT>) 

+ 4 ~ 3 (] tmu> -i mur> +1 rum> -i mm> +i urm> -i mru>)> i 

This state is odd under P-parity. The spectrum of the six sites chain is reported in fig.(^). There 
are 9 triplets in the spectrum. In J37J it was already pointed out that the number of lowest lying 
triplets for a finite system with N sites is N(N + 2)/8, so for N = 6 there are 6 lowest lying triplet 
states. In order to identify these 6 states among the 9 that are exhibited by the spectrum of fig.(H), 
it is necessary to compute their A's and their Q's. In this way in fact, we can find out which arc the 
triplets characterized by two holes and thus belonging to the triplet of type ( |S.158 ) . In table (|J) we 
report the internal quantum numbers of the lowest lying triplets. The Q^ h ''s vary in the segment 
[—§,§]■ The highest weight state of the triplet of zero momentum and energy — (J/2)(5 + v5) reads 



J 



0,0,--(5 + V5),0> 



1 



-A- 



-3 + V5 



\/45 - 15\/5 2 

i mm>+i tum>) 

(i mm> +i tmu> +i mm> +i mm> 
(i - V5)(i mm> +i mm> +i umt>)} 



(I Utttt> +1 mm> +1 TTTTii> +1 TTTiit> 



.196) 



tnnt> 



nntt>) 

(8.197) 



One can get the triplet of energy 
be explicitly checked from 



(J/2) (5 - \/5) from (8.197) by changing y/E -> ~VE. As can 



197), the two non-degenerate triplets of zero momentum are then 



P-parity even, namely they have opposite parity with respect to that of the ground state, as it 
happens for the lowest lying triplet excitations of the 2-flavor Schwinger model. For what concerns 
the degenerate triplets of momenta ir/3 and 57r/3 (or 27r/3 and Att/3) they do not have definite 
P-parity, but it is always possible to take a linear combination of them with parity opposite to the 
ground state. 



Table 3: Triplet internal quantum numbers 



TRIPLET Ai 



1 1 5+V5 t n y / 5-2VE 5-2V5 _1 1 

' '2 J '"^ V20 V 20 2 2 



1 1 5-V5 j n ^ 1 5+2V5 / 5+2 VZ _3 3 

x ' x 1 2 u i u ^ V20 V 20 2 2 



1 1 T TL ~> \pj+\fTT v3— \pii 3 1 

' ' 2 5 3 8 8 2 2 



1 1 _I+VH 7 2tt -> -2^-^-2+2717 -2V3+V-2+2yT7 _3 _1 
' ' 4 J ' 3 ^ 2+2VT7 2+2VT7 2 2 



1 1 T+VV7 j 4tt ^ 2x/3-V-2+2y / 17 273+^-2+2^ 1 3 

' ' 4 J ' 3 ^ 2+2VT7 2+2VTF 2 2 

1 1 _5 t 57T ^ V3 — Vtt V^+y^r _1 3 

' ' 2 5 3 8 8 22 



The remaining three triplets in fig. (j^) have no real A's and are characterized by a string of length 
2 and four holes for Q = —§,—5, 5, §, i.e. do not belong to the type (8.158). More precisely, two 
triplets have a string approximately of length 2, due to the finite size of the system, while the triplet 
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of momentum tt has a string exactly of length 2. In table 
of these triplets. 



we summarize the quantum numbers 



Table 4: Four holes triplet internal quantum numbers 



TRIPLET 




Ai 


A 2 


IM, 7 -fj, 


2-ir . 
3 ^ 


2 v / 3-'i\/ 2 + 2 v / T7 
-2+2V17 


2V3+JV / 2+2v / Tf 
-2+2\/T7 


|1,1,-J,7T> 




i 
2 


( 

2 




4ir . 
3 ^ 


2v / 3+i\/ 2+2\/l7 


2^-1-^/2+2717 




2-2VT7 


2-2v / 17 



In fig. (j^) it is shown that the spectrum exhibits five singlet states. The lowest lying state 
at momentum tt is the ground state. Then there are three excited singlets characterized by the 
configuration with two holes ( 8.157 ), i.e. they have one real A and a string of length almost 2. 
In table (JsJ) we summarize their quantum numbers. Among these singlets, those which are not 
degenerate, have P-parity equal to that of the ground state (odd) as it happens in the 2-flavor 
Schwinger model. The non-degenerate singlet in fact reads 



lo,o,-3j,o > = mm> +1 ntut> +1 mm> +i mm> +i utut> +i mm> 

V 12 

- 1 uun> -1 mm> -i tmu> -i mut> -i mm> -i tutit>} (&.198) 

The degenerate singlets are again not eigenstates of the P-parity, but it is always possible to take a 
linear combination of them with the a P-parity that coincides with that of the representative state 



198) of the configuration. 



Table 5: Singlet internal quantum numbers 



SINGLET 


A 


As 




|0,0,-3J,0 > 








-1 1 


|0,0,-2J,f > 


V3+2VE 

14 


-2+3V2 

V3(4+v^) 


1 


|0,0,-2J,f > 


\/3+2\/6 


2-3V2 


-1 


14 


V3(4+V2) 



The remaining singlet |0,0, — - 



— J, 7r > it is not of the type (8.157). It is characterized by a 

and A3,i = - 



- -•- /5+2V13 



12 



A 



2.1 



5+2V13 
12 



string approximately of length 3 with Ai 

Even in finite systems very small like the 4 and 6 sites chains, the "string hypothesis" is a very 
good approximation and it allows us to classify and distinguish among states with the same spin. 

The ground state of the antiferromagnetic Heisenberg chain with N sites is a linear combination 



N 

N 
2 



of all the ( Vr I states with spins up and spins down. These states group themselves into 



27 
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2tt/3 



7T 



4tt/3 



5tt/3 



2tt 



-1/2 

-(5- vT3)/2 

-(V- ViV)/2~ 



e 
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-5/2 



-e- 



-(7l7 + 7)/4 



-e- 



■o 



-0- 



■G 



-(V5 + 5)/2 



-(v / T3 + 5)/2 
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Figure 3: Six sites chain spectrum 



sets with the same coefficient in the linear combination according to the fact that the ground state 
is translationally invariant (with momentum (n) for ^ even (odd)), it is an eigenstate of P-parity 
and it is invariant under the exchange of up with down spins. The states belonging to the same set 
have the same number of domain walls, which ranges from N, for the two Neel states, to 2 for the 
states with ^ adjacent spins up and ^ adjacent spins down. 

The ground state of the 8 sites chain is 



\g.a. >= 
where 



8 > +a\4 1] > +/W > +7l< J > +<W > +^4 V > +CIV>2 >) (8.199) 



,,(2) 



(2) 



,,(3) 



iv* >=innnn>+i mum> 

I TTIT4TII> +1 44T4T4TT> ^translated states 

I TT4T44T4> +1 44T4TT4T> ^translated states 
I TT44TT44> ^translated states 

I TTT444T4> +1 Ulttt.Lt> ^translated states 

I TT444TT4> +1 UtttUt> ^translated states 
TTTT4444> ^translated states . 



>- 

\4 2) >- 
\*P >- 
\4 2) >- 
\4 3) >= 

|V>2 > = 



By direct diagonalization one gets 

a = -0.412773 



(8.200) 
(8.201) 
(8.202) 
(8.203) 
(8.204) 
(8.205) 
(8.206) 

(8.207) 
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(3 = 0.344301 

7 = 0.226109 

6 = -0.087227 

e = 0.136945 

C = 0.018754 

N = 2 + 16a 2 + 8/3 2 + 4 7 2 + 16<5 2 + 16e 2 + 8C 2 = 6.30356 



8.208) 
8.209) 
8.210) 
8.211) 
8.212) 
8.213) 



The energy of the ground state is 

Eg. s . = -5.65109J . (8.214) 



Eq. fl3.214 ) differs only by 1.8% from the thermodynamic limit expression E g . s . — — 81n2 = —5.54518 



Moreover also the correlation function of distance 2 Eq.(5.82) computed for the 8 sites chain is 
G(2) = 0.19577V, value which is 7% higher than the exact answer Eq. fl5.83p . 



In the analysis of finite size systems we were able to find the coefficient f3 of the first set of states 
containing N — 2 domain walls in the ground state. These states are obtained interchanging two 
adjacent spins in the Neel states. The f3 is for a generic chain of TV-sites 

N + 2E 0S , 
f3= g - s - = l-2m2 . (8.215) 

Appendix B: The correlator G(2) in terms of spin configura- 
tion probabilities 

In this appendix we shall establish a relation between the Heisenberg model correlator G(2) 
and the probabilities of finding, in the antiferromagnetic vacuum, certain groups of spin in a given 
position. 

The isotropy of the Heisenberg model implies that 

N N 

Y,<9-S-\S x -S x+2 \g. S .>=3Y / <9-S-\S 3 x -S 3 x+2 \g.s.> . (8.216) 

x — l x—l 

Let us introduce the probability P3 for finding three adjacent spins in a given position in the 
Heisenberg antiferromagnetic vacuum. Taking advantage of the isotropy of the Heisenberg model 



ground state and of its translational invariance, it is easy to see that the correlator (8.21£) can be 
written in terms of the P3S as 

N 1 

< 9-s.\S 3 x ■ S 3 x+2 \g.s >= N - 2 ( P 3 (TTT) + ft(UT) - ^(TTI) - JMlTt) ) ■ (8.217) 

x—l 



The factor 2 appears in (8.217) due again to the isotropy of the Heisenberg model: the probability 
of a configuration and of the configuration rotated by tt around the chain axis, are the same. 

In 111] the so called "emptiness formation probability" P(x) was introduced. 

X 

P(x)=<g.s.\l[P j \g.s.> , (8.218) 
3=1 

where 

Pj = 1(^ + 1) (8.219) 

and a 3 is the Pauli matrix. P{x) determines the probability of finding x adjacent spins up in the 
antiferromagnetic vacuum. One gets 

P(|||) = p(3) (8.220) 

P(T4T) = P(l) - 2P(2) + P(3) (8.221) 

P(4TT) = P(TT4) = P(2)-P(3) (8.222) 
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so that Eq.(5.82) reads 



^ = 2P(3) - 2P(2) + ~P(1) . (8.223) 



Using the exact value the correlator G(2) computed in [|17| from ( J8.223 ) and from the known values 
of P(2) and P(l) given in |l§ 

P(l) = 1 (8.224) 
P(2) = -(l-ln2) (8.225) 



.226) 



we get 



P(3) = i-ln2 + ^C(3) • (8.227) 
For the general emptiness formation probability P(x) of the antiferromagnetic Heisenberg chain, 



an integral representation was given in [IS]], but, to our knowledge, the exact value of P(3) (8.227) 
was not known. 
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